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We consider a question motivated by the third law of thermodynamics: can there be a local 
temperature arbitrarily close to absolute zero in a nonequilibrium quantum system? We consider 
nanoscale quantum conductors with the source reservoir held at finite temperature and the drain 
held at or near absolute zero, a problem outside the scope of linear response theory. We obtain 
local temperatures close to absolute zero when electrons originating from the finite temperature 
reservoir uirdergo destructive quantum interference. The local temperature is computed by numer¬ 
ically solving a nonlinear system of equations describing equilibration of a scanning thermoelectric 
probe with the system, and we obtain excellent agreement with analytic results derived using the 
Sommerfeld expansion. A local entropy for a nonequilibrium quantum system is introduced, and 
used as a metric quantifying the departure from local equilibrium. It is shown that the local entropy 
of the system tends to zero when the probe temperature tends to zero, consistent with the third law 
of thermodynamics. 
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I. INTRODUCTION 

The local temperature of a quantum system out of 
equilibrium is a concept of fundamental interest in 
nonequilibrium thermodynamics. Out of equilibrium, the 
temperatures of different degrees of freedom generally do 
not coincide, so that one must distinguish between mea¬ 
sures of lattice temperatureic— , photon temperature^^— , 
and electron temperature^r— . The Scanning Thermal 
Microscope (SThM)ii couples to all these degrees of free¬ 
dom, and thus measures some linear combination of their 
temperatures^^ in the linear response regime. Recent 
advances in thermal microscopyi^^— have dramatically 
increased the spatial and thermal resolution of SThM, 
pushing it close to the quantum regime. 

In this article, we focus on the local electron tempera¬ 
ture Tp as defined by a floating thermoelectric probei^^— . 
The probe, consisting of a macroscopic reservoir of elec¬ 
trons, is coupled locally and weakly via tunneling to the 
system of interest, and allowed to exchange charge and 
heat with the system until it reaches equilibrium, thus 
defining a simultaneous temperature and voltage mea¬ 
surement. Several variations on this measurement sce¬ 
nario have also been discussed in the literature^r— . 

The above definition of Tp is operational: tempera¬ 
ture is that which is measured by a suitably defined ther¬ 
mometer. Nonetheless, it has been shown that this defini¬ 
tion of Tp is consistent with the laws of thermodynamics 
under certain specified condition a^^d^ . In the present ar¬ 
ticle, we show that Tp is consistent with the third law of 
thermodynamics. In particular, we introduce a definition 
of the local entropy Ss of a nonequilibrium quantum sys¬ 
tem, and show that S's 0 as Tp ^ 0. Moreover, we 
show that values of Tp arbitrarily close to absolute zero 
can exist in quantum systems under thermal bias. Sg is 
also used to quantify the departure from local equilibrium 
beyond the linear response regime. 


The article is organized as follows: The nonequilibrium 
Green’s function (NEGF) formalism needed to describe 
the local properties of a nonequilibrium quantum system 
and its interaction with a scanning thermoelectric probe 
is introduced in Sec. m The local entropy of a nonequilib¬ 
rium quantum system is defined in Sec. Ill Bl along with a 
normalization that takes into account spatial variations 
in the local density of states. Analytical results for the 
minimum local temperatures in quantum systems under 
thermal bias are derived in Sec. nil The local temper¬ 
ature and entropy distributions in several 7r-conjugated 
molecular junctions under thermal bias are computed in 
Sec.[IVl Our conclusions are summarized in Sec. El while 
some useful details of the formalism and modeling are 
provided in Appendices lAl ICl 


II. FORMALISM 

We consider a temperature/voltage probe coupled lo¬ 
cally (e.g., via tunneling) to a nonequilibrium quantum 
system. The probe is also connected to an external 
macroscopic reservoir of noninteracting electrons held at 
a fixed temperature Tp and chemical potential ytp. We 
use the NEGF formalism to write the electron number 
current and heat current flowing into the probe as 

/ OO 

, 1 ) 

Tr {r>'(i.,)(G'(^) + - G‘‘(i.,)l)} . 

where v = Q refers to the electron number current— 
and V = \ gives the electronic contribution to the heat 
current—. G'^{u)) and G°‘{uS) are the Fourier transforms 
of the retarded and advanced Green’s functions, respec¬ 
tively, describing propagation of electronic excitations 
within the system, and G^{u}) is the Fourier transform 
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of the Keldysh “lesser” Green’s function describing the 
nonequilibrium population of the electronic spectrum of 
the system. rP(a;) is the tunneling-width matrix de¬ 
scribing the coupling of the probe to the system and 
fp{uj) = 1/(1 + exp )) is the equilibrium Fermi- 
Dirac distribution of the probe. Eq. o is a general re¬ 
sult valid for any interacting nanostructure under steady- 
state conditions. 


A. Local temperature and voltage measurements 

A dehnition for a local electron temperature and volt¬ 
age measurement on the system that takes into account 
the thermoelectric corrections was proposed in Ref. [l^ by 
noting that the temperature Tp and chemical potential ^p 
should be simultaneously defined by the requirement that 
both the electric current and the electronic heat current 
into the probe vanish: 


Eqs. allow us to rewrite Eq. o in a form analogous 
to the two-terminal Landauer-Biittiker formula 

27rTr{rP(a;)A(a;)} [fs{uj) - /p(w)]. 

It was noted that, for the case of maximum local cou¬ 
pling, [r^’(a;)]y = TP{uj)6inSjn, the quantities A{uj) = 
Ann(^) and fs{uj) become independent of the probe cou¬ 
pling, and can be related by the familiar equilibrium- 
type formula, = 27riA„„(a;)/s(a;), even though 

the system is out of equilibrium. 

In the broad-band limit rP(w) ~ r^(^o): where is 
the equilibrium Eermi energy of the system, we may write 
Tr{rP(^o)} = fP so that A{uj) = Tr {rP(/xo)A(a;)}/fP. 
Erom Eq. 0, we have 

pOO 

Ip '"''=Yy ~ (8) 


=0, 1/ e {0,1}. (2) 

Eq. (21) gives the conditions under which the probe is in 
local equilibrium with the sample, which is itself arbi¬ 
trarily far from equilibrium. 

Previous analysesi^^— have considered this problem 
within linear response theory, which reduces the system 
of nonlinear equations ([2]) to a system of equations linear 
in Tp and ^p. In this article, we consider a problem that 
is essentially outside the linear response regime and solve 
the nonlinear system of equations m numerically. 

It was shown in Ref. that Eq. o can be written 
in terms of the local properties of the nonequilibrium 
system. The mean local spectrum sampled by the probe 
was defined as 


A{uj) 


Tr{rP(a;)A(w)} 

Tr{rP(w)} 


( 3 ) 


where A{uj) = i{G^{uj) — G“(w))/27r is the spectral func¬ 
tion of the nonequilibrium system. Motivated by the re¬ 
lation at equilibrium, = 2TriA{uj)feq{uj), the lo¬ 

cal nonequilibrium distribution function (sampled by the 
probe) was defined as 


^ Tr{rP(a;)G<(a;)} 
27ri Tr {rP(a;)A(w)} 


( 4 ) 


The mean local occupancy of the system orbitals sampled 
by the probe is^^ 



duA{uj)fs{ui), 


( 5 ) 


and similarly, the mean local energy of the system or¬ 
bitals sampled by the probe is^^ 



diOLU A{lij) fs{<-o). 


( 6 ) 


It was noted that the equilibrium condition of Eq. now 
implies that the mean local occupancy and energy of the 
nonequilibrium system are the same as if its nonequilib¬ 
rium spectrum A{u}) were populated by the equilibrium 
Fermi-Dirac distribution of the probe fp{uj): 

{»>!,, = («) I,. (9) 

= WI/.. (“) 

i.e., the probe equilibrates with the system in such a way 
that fp{uj) satisfies the constraints imposed by Eqs. dU 
and (fTUl) . 


B. Local entropy 

The von Neumann entropy^^ for a system of noninter¬ 
acting fermions can be expressed in terms of the single¬ 
particle occupation probabilities pi as 

S = -Y^[pAnp^ + (1 -p01n(l -p^)], (11) 

i 

which can be extended to the case of a continuous spec¬ 
trum by simply replacing the summation by an integral 
over the density of states. We propose a natural exten¬ 
sion for the “local entropy” Ss of the nonequilibrium sys¬ 
tem, within an effective one-body description, with the 
local density of states (sampled by the probe) given by 
A{uj), and the occupation probabilities given by the local 
nonequilibrium distribution of the system fs- 

/ OO 

dujA{uj)[fs{oj) lnfs{oj) 

-OO 

+ (1 - /s(w))ln(l - /s(w))]. 

( 12 ) 

Within elastic transport theory, it can be shown that 
0 < /s < 1 1 Appendix lAT) and therefore Ss in Eq. (fT^ 
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is real and positive. Ss also correctly reproduces two 
known limiting cases for the entropy of a system of inde¬ 
pendent fermions: (i) in equilibrium, gives the correct 
entropy of the subsystem sampled by the probe; (ii) Sg 
gives the correct nonequilibrium entropy for an entire sys¬ 
tem of fermions^. However, we note that the entropies 
(na of the various subsystems of a quantum system are 
not additive out of equilibrium. The definition of local 
nonequilibrium entro py given by Eq. (El) differs from 
that proposed in Ref. [^. 

We also define the local entropy of the corresponding 
local equilibrium state of the system if its local spectrum 
were populated by the probe’s equilibrium distribution 
function fpi 

/ OO 

du}A{u})[fp{uj)\n fp{u}) 

-OO 

+ (1 - /p(w))ln(l - /p(w))]. 

(13) 

For sufficiently low probe temperatures, 

Sp ~ ^A{fio)kBTp, (14) 

a standard textbook result. The maximum entropy prin¬ 
ciple implies Sp > Sg, since the Fermi-Dirac distribu¬ 
tion /(w) = fp{uj) maximizes the local entropy S'[/(u;)] 
subject to the constraints imposed by Eqs. dH) and (El. 
Clearly, 5^ —>■ 0 as Tp —>■ 0, which implies Sg ^ 0 as 
Tp —>■ 0 (third law of thermodynamics). 

We propose the local entropy deficit AS = Sp — Sg 
as a suitable metric quantifying the departure from lo¬ 
cal equilibrium. However, it is important to note that 
the mean local spectrum A(uj) varies significantly from 
point to point within the nanostructure depending upon 
the local probe-system coupling (especially in the tun¬ 
neling regime) and limits the use of AS* while compar¬ 
ing the ‘distance’ from equilibrium for points within the 
nanostructure. The situation is analogous to that of 
a dilute gas, which can have a very low entropy per 
unit volume even if it has a very high entropy per par¬ 
ticle. We note that states far from the equilibrium 
Fermi energy /Tp contribute negligibly to the entropy since 
lim/_j.o-S')/] = lim/^i S[f] = 0, and therefore introduce 
a normalization averaged over the thermal window of the 
probe: 


Ar = 


duj 


A{uj) f-df, 


Tr{H(w)}V duj 


(15) 


We define the local entropy-per-state of the system Sg and 
that of the corresponding local equilibrium distribution 


Sr, as 


Sg 

Sp 

^^ = 77- 


(16) 

(17) 


As = Sp — Sg quantifies the per-state ‘distance’ from lo¬ 
cal equilibrium. We present numerical calculations of the 
local entropy-per-state in Sec. IIVI and discuss its impli¬ 
cations. 


III. LOCAL TEMPERATURES NEAR 
ABSOLUTE ZERO 

Our analyses in this article consider a quantum con¬ 
ductor that is placed in contact with two electron reser¬ 
voirs: a cold reservoir R1 and a hot reservoir R2. We 
are interested in the limiting case where reservoir R1 is 
held near absolute zero (Ti —> 0) while R2 is held at fi¬ 
nite temperature (T 2 = 100A in our simulations). We 
assume no electrical bias (/ii = p ,2 = Mo)- 

Transport in this regime will be dominated by elastic 
processes, and occurs within a narrow thermal window 
close to the Fermi energy of the reservoirs. It should 
be noted that, although the transport energy window is 
small, the problem is essentially outside the scope of lin¬ 
ear response theory owing to the large discrepancies in 
the derivatives of the Fermi functions of the two reser¬ 
voirs. Therefore, the problem has to be addressed with 
the full numerical evaluation of the currents given by Eq. 

(ED- 

For many cases of interest, transport in a nanostruc¬ 
ture is largely dominated by elastic processes. This al¬ 
lows us to utilize the simpler formula analogous to the 
multiterminal Biittiker formula^ given by— 

1 _ foo 

duj{uj - iip)''TpA^){fa{uj) - fv{^)), (18) 

n ^ 7-00 

where the transmission function is given by^ 

Tpa^iu:) = Tr{rP(w)G'-(a;)r“(w)G“(w)}. (19) 

A pure thermal bias, such as the one considered 
in this article, has been shown to lead to temper¬ 
ature oscillations in small molecular junctionsi^ and 
1-D conductor s'^ 1 ^^ . Temperature oscillations have 
been predicted in quantum coherent conductors such as 
graphenei^, which allow the oscillations to be tuned (e.g., 
by suitable gating) such that they can be resolved un¬ 
der existing experimental techniques'^— of SThM. More 
generally, quantum coherent temperature oscillations can 
be obtained for quantum systems driven out of equi¬ 
librium due to external fields^i^ as well as chemical 
potential^! and temperature bias of the reservoirs. In 
practice, the thermal coupling of the probe to the envi¬ 
ronment sets limitations on the resolution of a scanning 
thermoelectric probe^^. However, in this article, we ig¬ 
nore the coupling of the probe to the ambient environ¬ 
ment, in order to highlight the theoretical limitations on 
temperature measurements near absolute zero. 

In evaluating the expressions for the currents in Eq. 
(1181) within elastic transport theory, we encounter inte¬ 
grals of the form j7°^dujF{u})(^f2{uj) — fi{ui)). We use 
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the Sommerfeld series given by 


Ti —> 0 which yields 


/ oo 

dujF{u!)[f2{u!) — = / dijF{uj) 

- OO J 

+ 2 ^ e(fc + 1 ) [{kBTif+^F^’^'' {^, 2 ) - (fesTO'^+'F''') ifii)], 

k 

ke {1,3,5,...}, 

( 20 ) 

where we use the symbol 0 that relates to the Riemann- 
Zeta function as0(fc+l) = (l — ^)C(fc + l) and /^(a;) is 
the Fermi-Dirac distribution of reservoir a. The second 
term on the r.h.s of Eq. (1201) accounts for the exponential 
tails in (/ 2 (w) — /i(w)), and its contribution depends on 
the changes to the function F{ui) in the neighbourhoods 
of oj = /ii and oj = ^2 and can generally be truncated 
using a Taylor series expansion for most well-behaved 
functions F{oj). The l.h.s of Eq. (1^ is bounded if F{uj) 
grows slower than exponentially for uj —>■ ±00 and is sat¬ 
isfied by the current integrals in Eq. dUD. 


A. Constant trasmissions 


T = 


% 


p2 


7pi + Tp2 


To. 


(26) 


We have Tp —)> 0 as Tp 2 0. Indeed, when the system 
is decoupled from the hot reservoir R2, the probe would 
read the temperature of reservoir Rl. 


B. Transmission node 

The analysis of the previous section suggests that 
a suppression in the transmission from the finite- 
temperature reservoir R2 results in probe temperatures 
in the vicinity of absolute zero. In quantum coherent 
conductors, destructive interference gives rise to nodes 
in the transmission function. In this section, we consider 
a case where the transmission from R2 into the probe 
has a node at the Fermi energy. In the vicinity of such 
a node, generically, the transmission probability varies 
quadratically with energy: 


In order to make progress analytically, we consider first 
the case of constant transmissions: 


= (27) 


Fpcxij-O') — 'Tpa^Fo') — 'Tpa- 


( 21 ) 


while the transmission from the cold reservoir i?l may 
still be treated as a constant: 


This is a reasonable assumption because the energy win¬ 
dow involved in transport is of the order of the thermal 
energy of the hot reservoir {kBT 2 ~ 25meV, at room 
temperature) and we may expect no great changes to the 
transmission function. In this case the series (|20)) for the 
number current contains no temperature terms at all, 
while the heat current contains terms quadratic in the 
temperature. It is easy to see that the expression for the 
number current into the probe becomes 

4°^ = ( 22 ) 

Q: 

and the heat current into the probe is given by 


Tpi{oj) — Tpi. (28) 

Applying the Sommerfeld series (I20p for the number cur¬ 
rent gives us 


1 / 

4°^ =T - Mp) Mof 

n\ b ^ 29 ) 


where the fc_BT 2 term is still missing since the first deriva¬ 
tive of Tp 2 {oj) vanishes at ^2 = Mo- We note that Eq. (|2^ 
admits a single real root at 


Q, \ / 

(23) 

Eq. (j22l) does not depend on the temperature and can be 
solved readily: 

dp = do, (24) 

since di = d 2 = do and Eq. (USD is solved by 


Tp = 


TplTi + Tp2T2 
Tpi + Fp2 


(25) 


In this article, we are primarily interested in temperature 
measurements near absolute zero and work in the limit 


flp — ^ 0 - 


(30) 


With this solution, we can write down the equation for 
the heat current as 


7 ( 1 ) 

p 
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8 15 



which gives us a simple quadratic equation in T^. We 
note that the above equation is monotonically decreasing 
in Tp for all positive values of temperature. There exists 
a unique solution to Eq. m in the interval Ti < Tp < 
T 2 , since 4^\^p) undergoes a sign change between these 
two values, and is also the only positive solution due 
to monotonicity. Physically, this solution is reasonable 
since we expect a temperature measurement to be within 
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the interval (Ti,T 2 ) in the absence of an electrical bias. 
It is straight-forward to write down the exact solution 
to Eq. (IMT) (see Appendix ICl) . However, we simplify the 
expression for Tp by noting that 


T^?{kBT2f <^Tpu (32) 


that is, the transmission into the probe from R2 within a 
thermal energy window kBT 2 in the presence of a node, 
is small in comparison to the transmission from i?l. The 
approximate solution for Tp then becomes 


Tp = 


h^^ri;^{kBT2)\ 


20 


Tpi 


(33) 


where, as before, we have taken the limiting case where 
Ti 0. 


C. Higher-order destructive interferences 

Although a generic node obtained in quantum coher¬ 
ent transport depends quadratically on the energy, it 
is possible to obtain higher-order “supernodes” in some 
systems^. In the vicinity of such a supernode, the trans¬ 
mission function can be written as 

= (.-«)’”, (34) 

while the transmission from i?l may still be approxi¬ 
mated by Eq. (l28l) . Exact expressions for the currents 
can again be evaluated using Eq. (EOl). The expression 
for the number current becomes 

. X .T-( 2 n) 

+2 ^ e(fc+l)[(fcBT 2 )"+V«(ho)-(fcnrp)'=+'r«(Mp)]), 

k^odd 

(35) 

and we have 

r^2\t^o)=0, VfcG{l,3,5,...}. (36) 

Now, with 

k-p ~ kOi ( 3 "^) 

every term on the r.h.s of Eq. (ESI) vanishes. With this 
solution for /ip, we proceed to write the equation for the 
heat current. Using Eq. ( 1201 ) with F{uj) = {uj - fJ.o)Tpa{uj) , 
we obtain only one nonvanishing derivative for each reser¬ 
voir, that is, U^^^(/io) = Tpiifio) for i?l and E^^"+^^(/io) = 
(2n -I- 1 ) 7 ^ 2 ^"^ for R2. Therefore, 

4 '’ =|(e(2)7;i[(fcBTi)" - [ksTpf] 

-t (2n + l)e(2n -t [{kBT 2 f"+^ - {kBTpf^+^]J, 

(38) 


which is a polynomial equation in Tp of degree ( 2 n-|- 2 ). 
We can rewrite Eq. (1551) as a polynomial p{x) in x = 
TplT2: 

p{x) = x'^ + - A„ = 0, (39) 


where we have taken Ti —^ 0 for i?l, and A„ is a dimen¬ 
sionless quantity given by 


( 2 n+l) 


e{2n + 2) 
0 ( 2 ) 




(40) 


We will have A„ <C 1 for a suitable energy window set 
by fc_BT 2 , since the transmission into the probe from R2 
suffers destructive interference at the Eermi energy (/xg). 
If the thermal energy is large enough, then this approx¬ 
imation may no longer hold. In any case, it is possible 
to define a temperature T 2 so that this approximation is 
strongly valid. Under the validity of this approximation, 
the solution to Eq. (1391) can be written using perturbation 
theory as 

x=./k(^1 + 0{X^+^)'^, (41) 

with corrections that are of much higher-order in A„. The 
solution for Tp given in Eq. (1411) reduces to Eq. (|26|) for 
the case with constant transmissions by setting n = 0 , 
and to the approximate result obtained in Eq. (1331) in the 
presence of a node by setting n = 1. We note that higher- 
order interference effects cause the probe temperature to 
decay more rapidly with respect to T 2 since Tp 
that is, when T 2 is halted, Tp is reduced by a factor of 
2"+^. Now, if we consider the limiting case where R2 is 
also cooled to absolute zero, T 2 —?► 0, Eq. (ITT]) implies 
that Tp ^ 0 at least as quickly as T 2 (in the absence 
of destructive interference, i.e., n = 0 ) or quicker (when 
there is destructive interference, i.e., n > 1 ). 

It should be noted that the polynomial given in Eq. 
( 1391 ) is monotonically increasing for all positive x and fur¬ 
thermore, there is only one positive root since p( 0 ) < 0 
and p(l) > 0. Stated in terms of Tp, this implies the exis¬ 
tence of a unique solution for the measured temperature 
Tp in the interval Ti < Tp < T 2 , as noted in the previous 
subsection (Ti has been set to zero in Eq. ( 1551 ) 1 . The sign 
of p{x) essentially tells us the direction of heat flow for 
a temperature bias of the probe with respect to its equi¬ 
librium value. p{x) < 0 (j>{x) > 0) corresponds to heat 
flowing into (out of) the probe, since we changed the sign 
in re-arranging Eq. ( 1551 ) to Eq. ( 1591 ) . The monotonicity 
of p{x) is therefore equivalent to the Clausius statement 
of the second law of thermodynamics, and also ensures 
the uniqueness of the temperature measurement. 


IV. NUMERICAL RESULTS AND DISCUSSION 

In order to illustrate the above theoretical results, and 
further characterize the local properties of the nonequi¬ 
librium steady state, we now present numerical calcula¬ 
tions for several molecular junctions with tt- conjugation. 
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In all of the simulations, the molecule is connected to a 
cold reservoir R1 at Ti = OK and a hot reservoir R2 at 
T 2 = lOOK. There is no electrical bias; both electrodes 
have chemical potential /tq. The temperature probe is 
modeled as an atomically-sharp Au tip scanned horizon¬ 
tally at a constant vertical height of S.SA above the plane 
of the carbon nuclei in the molecule (tunneling regime). 

The molecular Hamiltonian is described within Hiickel 
theory, i/inoi = tydjdj + h.c, with nearest-neighbor 
<i,j> 

hopping matrix element t = —2.7eV. The coupling of 
the molecule with the reservoirs is described by the 
tunneling-width matrices r“. The retarded Green’s func¬ 
tion of the junction is given by G'^{uj) = [Stu — i/moi — 
where tunneling 

self-energy. We take the lead-molecule couplings in the 
broad-band limit, i.e., r“^(a;) = r“^(//o) where /xq is the 
Fermi energy of the metal leads. We also take the lead- 
molecule couplings to be diagonal matrices r“^(w) = 
^a^ni^mi coupled to a single 7r-orbital I of the molecule. 
S is the overlap-matrix between the atomic orbitals on 
different sites and we take S = I, i.e., an orthonormal 
set of atomic orbitals. The lead-molecule couplings are 
taken to be symmetric, with Ti = r 2 = 0.5eV. The non¬ 
zero elements of the system-reservoir couplings for i?l 
(cold) and R2 (hot) are indicated with a blue and red 
square, respectively, corresponding to the carbon atoms 
in the molecule covalently bonded to the reservoirs. The 
tunneling-width matrix T^* describing probe-sample cou¬ 
pling is also treated in the broad-band limit (see Ap¬ 
pendix m for details of the modeling of probe-system 
coupling). The probe is in the tunneling regime and the 
probe-system coupling is weak (few meV) in comparison 
to the system-reservoir couplings (Ti = r 2 = O.SeV). 

It must be emphasized that, although we take a non¬ 
interacting Hamiltonian for the isolated molecule, our 
results depend only upon the existence of transmission 
nodes, which are a characteristic feature of coherent 
transport, and do not depend on the particular form of 
the junction Hamiltonian. 


A. Local temperatures 

We considered several different molecules and electrode 
configurations, with and without transmission nodes, and 
searched for the coldest spot in each system (indicated by 
a green dot in the figures) as measured by the scanning 
thermoelectric probe. 

The local temperature depends on the transmissions 
from the reservoirs into the probe (Eq. (1191) 1. determined 
by the local probe-system coupling T^* (see appendix IbI) . 
and is thus a function of probe position. The coldest 
spot was found using a particle swarm optimization tech¬ 
nique that minimizes the ratio of the transmissions to 
the probe (within a thermal window) from the hot reser¬ 
voir R2 to that of the cold reservoir Rl, within a search 
space that spans the z-plane at 3.5A and restricted in the 


xy direction within lA from the edge of the molecule^. 
The numerical solution to Eq. @ was found using New¬ 
ton’s method. While the algorithm was found to con¬ 
verge rapidly for most points (less than 15 iterations), 
it is still computationally intensive since the evaluation 
of the currents given by Eq. dn must have sufficient 
numerical accuracy. We also note that the minimum 
probe temperature obtained for each junction does not 
depend strongly on the distance between the plane of 
the scanning probe and that of the molecule. This is ex¬ 
plained as follows: the probe temperature must depend 
upon the relative magnitudes of the transmissions into 
the probe from the two reservoirs, and not their actual 
values. Therefore, the temperature remains roughly in¬ 
dependent of the coupling strength TrjFP}. It has also 
been previously noted in Ref. El that the local tempera¬ 
ture measurement showed little change with the coupling 
strength even when varied over several orders of magni¬ 
tude. The restrictions placed on our search space within 
the optimization algorithm are therefore well justified. 

Fig. m shows the temperature distribution for two con¬ 
figurations of Au-benzene-Au junctions with the chem¬ 
ical potentials of the metal leads at the middle of the 
HOMO-LUMO gap. The mid-gap region is advanta¬ 
geous, since (i) the molecule is charge neutral when the 
lead chemical potentials are tuned to the mid-gap energy, 
and (ii) the mismatch between the metal leads’ Fermi en¬ 
ergy and the mid-gap energy is typically small (less than 
1-2 eV for most metal-molecule junctions) and available 
gating techniques^ would be sufficient to tune across the 
gap. In both junctions, the region of lowest temperature 
passes through the two sites in meta orientation relative 
to the hot electrode, because the transmission probabil¬ 
ity from the hot electrode into the probe is minimum 
when it is at these locationsi^. The meta junction ex¬ 
hibits much lower minimum temperature measurements 
since there is a transmission node from the hot reservoir 
R2 at the mid-gap energy, but such nodes are absent in 
the para junction. Table H] shows the coldest temperature 
found in each of the junctions presented here. 

We only present two junction geometries for benzene 
here to illustrate that the existence of a node in the probe 
transmission spectrum, at the mid-gap energy, depends 
upon the junction geometry. The nodes are also absent in 
the ortho configuration of benzene, and the lowest tem¬ 
perature found in this case is similar to that found in the 
para configuration. 

Fig. [5] shows the temperature distribution in a gated 
Au-pyrene-Au junction, and the transmissions into the 
probe from the two reservoirs at the coldest spot. In gen¬ 
eral, nodes in the transmission spectrum occur only in a 
few of the possible junction geometries. As in the benzene 
junctions, the coldest regions in the pyrene junction pass 
through the sites to which electron transfer from the hot 
electrode is blocked by the rules of covalenc^ describing 
bonding in 7r-conjugated systems. We note that the tem¬ 
perature distribution shown in Fig. [2] differs significantly 
from that shown in Ref. [l^ for four important reasons: (i) 
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FIG. 1: Upper panels: Local temperature distributions for Au-benzene-Au junctions in meta and para 
configurations, respectively. The thermal bias is supplied by cold (Ti = OK) and hot (T 2 = lOOK) reservoirs 
covalently bonded to the atoms indicated by the blue and red squares, respectively, and there is no electrical bias. 
The probe is scanned at a height of S.SA above the plane of the carbon nuclei in the molecule. The green dots 
shown in the temperature distributions correspond to the coldest temperature found in each of the junction 
configurations. Bottom panels: Transmission probabilities into the probe from i?l (cold, i.e., blue curve) and R2 
(hot, i.e., red curve), when the probe is positioned over the coldest spot (shown by the corresponding green dot in 
the upper panel). The existence of a transmission node in the meta configuration leads to a greatly suppressed 
probe temperature (see Table U for comparison). Note the very different vertical scales in the bottom panels. 


the junction configuration in Fig. [2] is asymmetric, while 
that considered in Ref. [I^ was symmetric; (ii) the ther¬ 
mal coupling Kpo of the temperature probe to the ambient 
environment has been set to zero in Fig.[2]to allow for res¬ 
olution of temperatures very close to absolute zero, while 
the probe in Ref. [l^ was taken to have Kpo = 
where Kpo = /3){kgT/h) = 2.84 x 10“^°W/K at 

T = 300K is the thermal conductance quantu m^^d^ ; (hi) 
the transport in Fig. [5] is assumed to take place at the 
mid-gap energy due to appropriate gating of the junction, 
while Re f. Il3l considered a junction without gating; and 
(iv) Ref. [l3| considered temperature measurements only 
in the linear response regime, while the thermal bias ap¬ 
plied in Fig. [2] is essentially outside the scope of linear 
response. Points (ii)-(iv) also differentiate the results 
for benzene junctions shown in Fig. [1] from the linear- 
response results of Ref. [H. 

Fig. [3] shows the temperature distribution in a gated 
Au-coronene-Au junction exhibiting a node in the probe 


transmission spectrum. The junction shown was one of 
three such geometries to exhibit nodes (10 distinct junc¬ 
tion geometries were considered). Again, the coldest re¬ 
gions in the junction pass through the sites to which 
electron transfer from the hot electrode is blocked by 
the rules of covalence. The coronene junction in Fig. [3] 
displays the lowest temperature amongst all the differ¬ 
ent junctions considered, with a minimum temperature 
of Tp = 35mK. It should be noted that this temperature 
would be suppressed by a factor of 100, i.e., Tp = 350^K 
if R2 were held at lOK due to the quadratic scaling of Tp 
with respect to the temperature T 2 [cf. Eq. (1331) ]. Higher- 
order nodes would produce even greater suppression. 

B. Local entropies 

The local temperature distributions shown in Figs. 
mm are essentially outside the scope of linear response 



































FIG. 2: Left panel: Probe temperature distribution in a Au-pyrene-Au junction under the same conditions described 
in Fig. [TJ The green dot corresponds to the coldest temperature found by the search algorithm. Right panel: 
Transmissions into the probe from the hot reservoir R2 (red) and the cold reservoir R1 (blue) at the coldest 
position, indicated by the green dot on the left. The probe transmission from R2 exhibits a (mid-gap) node at the 
Fermi energy fiQ of the reservoirs, thereby suppressing the temperature measured by the probe. 


Junction 

Tp[K] 

Numerical 

Tp[K] 

Analytic 

Eq. 

benzene (meta) 

0.154 

0.1526 

(|33|l 

benzene (para) 

4.624 

4.627 

(|26[ 

pyrene 

0.0821 

0.0817 

(|33|l 

coronene 

0.0349 

0.0355 

(|33|l 


TABLE I: The table shows the lowest temperatures 
found in the different junctions considered. All 
junctions have the same bias conditions: Ti = OK, 

T 2 = 100K and no electrical bias. The right-most 
column shows the equation used to compute the 
temperature analytically. We obtain excellent 
agreement between the numerical and analytic results. 
The para configuration of the benzene junction does not 
display a node in the probe transmission spectrum and 
therefore the minimum probe temperature is not 
strongly suppressed. 


theoryi^ since the cold reservoir R1 is held at Ti = OK, 
and derivatives of the Fermi function are singular at 
T = 0. However, it is an open question how far out of 
equilibrium these systems are and which regions therein 
manifest the most fundamentally nonequilibrium char¬ 
acter. To address such questions quantitatively, we use 
the concept of local entropy per state introduced in Sec. 
Ill Bl In particular, the normalized local entropy deficit 
As = Sp — Ss defined through Eqs. (inii-dni) allows us to 
quantify how far the system is from local equilibrium. 

Fig. S] shows the local entropy distribution of the sys¬ 
tem Ss and that of the corresponding local equilibrium 
distribution Sp, defined by Eqs. (USD and (CZD, respec¬ 


tively, for the Au-pyrene-Au junction considered earlier 
in Fig. [21 The Sp distribution strongly resembles the 
temperature distribution shown in Fig.|2l consistent with 
the fact ([T^ that the equilibrium entropy of a system of 
fermions is proportional to temperature at low tempera¬ 
tures. This resemblence is only manifest in the properly- 
normalized entropy per state Sp] the spatial variations 
of Sp are much larger, and stem from the orders-of- 
magnitude variations of the local density of states A{^o). 
The nonequilibrium entropy distribution Sg of the sys¬ 
tem qualitatively resembles Sp, but everywhere satisfies 
the inequality Ss < Sp (see Sec. IIIBIl . —S' 0 whenever 
Tp —>■ 0, consistent with the third law of thermodynam¬ 
ics. 

The deviation from local equilibrium is quantified by 
the local entropy deficit As = Sp — Ss shown in the top 
right panel of Fig.lD As shows deep blue regions (low en¬ 
tropy deficit) in both the hottest and coldest parts of the 
system, while the largest entropy deficits (bright red) oc¬ 
cur in the areas at intermediate temperatures. This may 
be explained as follows: within elastic transport theory, 
the local nonequilibrium distribution function is a linear 
combination of the distribution functions of the various 
reservoirs (see Appendix]^. The entropy deficit is min¬ 
imal when this distribution function strongly resembles 
the equilibrium Fermi-Dirac distribution of one of these 
reservoirs. Conversely, the entropy dehcit is maximal 
when there is a large admixture of both hot and cold elec¬ 
trons without inelastic processes leading to equilibration. 
Therefore, the hottest and coldest spots show the small¬ 
est entropy deficits since there is very little mixing from 
the cold reservoir R1 and hot reservoir R2, respectively, 
while the regions at intermediate temperatures have the 
largest entropy dehcitis, and hence are farthest from local 
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FIG. 3: Upper panels; Probe temperature distributions for a Au-coronene-Au junction under the same conditions 
described in Fig. [T] The numerically calculated temperature is on the left, and to the right is the analytically 
calculated temperature using Eq. (I26D . Although it is in excellent qualitative agreement (and quantitative 
agreement for the most part), Eq. (I26I) poorly estimates the temperature for the coldest spot (shown in green) due 
to the existence of a transmission node. Eq. (I33I) gives the correct estimation in the presence of a node, while Eq. 
(1^51) incorrectly predicts Tp = OK. The lower panel shows the probe transmissions from the two reservoirs (i?l with 
Ti = OK in blue and R2 with T 2 = lOOAT in red) corresponding to the probe positioned over the coldest spot (shown 
in green). 


equilibrium. 


V. CONCLUSIONS 


However, it can be seen that the colder spots are more 
strongly affected due to the mixing from the hot reser¬ 
voir R2, while the hotter spots are affected to a lesser 
extent due to the mixing from the cold reservoir Rl. 
This reflects the fact that the distribution function fs 
deviates much more from the distribution function /i (w) 
with Ti —> 0 (implying a pure state with zero entropy) 
due to a small admixture of hot electrons from R2 than 
is the case for the opposite scenario. In other words, it 
is easier to increase the entropy deficit As of a cold spot 
by adding hot electons (and thus driving it out of equi¬ 
librium) than it is the other way round. The entropy 
deficit is a good metric to capture such a change in the 
distribution function, and gives us a per-state ‘distance’ 
from local equilibrium. 


We investigated local electronic temperature distribu¬ 
tions in nanoscale quantum conductors with one of the 
reservoirs held at finite temperature and the other held 
at or near absolute zero, a problem essentially outside the 
scope of linear response theory. The local temperature 
was defined as that measured by a floating thermoelec¬ 
tric probe, Eq. ([2]). In particular, we addressed a question 
motivated by the third law of thermodynamics: can there 
be a local temperature arbitrarily close to absolute zero 
in a nonequilibrium quantum system? 

We obtained local temperatures close to absolute zero 
when electrons originating from the finite temperature 
reservoir undergo destructive quantum interference. The 
local temperature was computed by numerically solving 
a nonlinear system of equations [Eqs. ((2) and (flSl) ] de- 
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FIG. 4: Upper panels (left to right): The local entropy-per-state of the system Ss, of the corresponding local 
equilibrium distribution Sp, and the local entropy deficit As = Sp — Ss- The temperature distribution for the same 
junction (with identical bias conditions and sampling of probe positions) is shown in Fig. [21 and we note that it 
resembles almost exactly the distribution Sp. Lower panels: The distributions fs and fp for three points shown in 
the upper panels, each having different probe temperatures Tp = 0.15K, 15.9K and 81.6K, respectively. The 
corresponding entropy deficits are As = 0.11, 7.3, and 1.5, respectively, xl0“^. Point 2, although closer to OK than 
point 3 is to lOOK, is further from local equilibrium. 


scribing equilibration of a scanning thermoelectric probe 
with the system, and we obtain excellent agreement with 
analytic results [Eqs. (1^ . (1551) . and (HTl) ] derived using 
the Sommerfeld expansion. Our conclusion is that a lo¬ 
cal temperature equal to absolute zero is impossible in 
a nonequilibrium quantum system, but arbitrarily low 
finite values are possible. 
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A definition for the local entropy [Eq. (115]) ] of a 
nonequilibrium system of independent fermions was pro¬ 
posed, along with a normalization factor [Eq. (TTSl) ] that 
takes into account local variations in the density of states. 
The local nonequilibrium entropy is always less than or 
equal to that of a local equilibrium distribution with the 
same mean energy and occupancy, and the local entropy 
deficit was used to quantify the distance from local equi¬ 
librium in a nanoscale junction with nonlinear thermal 
bias (Fig. |4|). It was shown that the local entropy of the 
system tends to zero when the probe temperature tends 
to zero, implying that the local temperature so defined 
is consistent with the third law of thermodynamics. 


Appendix A: Elastic transport regime 

We derive the form of the nonequilibrium distribution 
function /s(a;) when the transport is dominated by elastic 
processes. We assume a nanostructure connected to M 
reservoirs, including the probe. Eq. o takes the form of 
Eq. (fT51) when the transport is elastic, and we have 

2^Tr{rP(u;)A(tc)}(/,(tc) - /^(u;)) 

M 

= Tpa{uj){fa{uj) - fp{uj)). 

Oi — 1 

(Al) 
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Now, we wish to rewrite the above equation in terms of 
the local properties sampled by the probe: 


Tr{rna.)A(c^)} 

Tr{rP(w)} 


- fp{u})) 


M 


= E 


Tr{rP(w)G’'(w)r“(w)G“(w)} 

27rTr{rP(w)} 


(/a(w) - /p(w)), 

(A2) 


where the Hrst factor on the l.h.s is the mean local spec¬ 
trum A{uj) sampled by the probe, defined by Eq. ([3|), 
and we used Eq. (IT^ for the elastic transmissions on the 
r.h.s. We dehne the injectivity of a reservoir a sampled 
by the probe as 


^ ^ ^ 1 Tr{rP(a;)G’'(w)r“(w)G“(a;)} 

PpaM - 2^ Tr{rP(w)} 


(A3) 


for the factors appearing on the r.h.s of Eq. (IA2() . In¬ 
jectivity of a reservoir a has been previously defined^ 
as the local partial density of states (LPDOS) associated 
with the electrons originating from reservoir a and, due 
to number conservation, the sum of injectivities of the 
reservoirs gives the local density of states (LDOS). We 
state an equivalent result for the injectivities defined in 
Eq. (|A3p in the following paragraph. Before proceeding, 
we note that the injectivities sampled by the probe, in 
Eq. (IA3I) . reduces to the LPDOS for electrons injected 
by reservoir a when the probe coupling is maximally lo¬ 
cal, i.e., [rP(w)]jj = TP{oj)Sin6jn and becomes essentially 
independent of the probe coupling when it is weak. Eq. 
(IA3I) also extends to a = p and defines the probe injec¬ 
tivity sampled by itself, which becomes negligible in the 
limit of weak coupling. 

It can be shown that the spectrum can be written as^ 


A{uj) = ^G’-(u;)r(cc)G“(a;), (A4) 

where r(a;) is given by 

r(a;)=^r“(w). (A5) 

Oi 

The contribution due to interactions r'"*(a;) in Eq. (IA5I) 
is missing since the interaction self-energy is Hermitian 
for elastic processes. Eqs. (IA3I) . (IA4I) and (IA5I) imply: 


M 

'^Ppa{i^) = A{uj). (A6) 

From Eq. (|A2p . we write 


M 

A{uj){fs{w)-fp{w)) = ^pp«(w)(/a(a;)-/p(w)) (A7) 
and Eq. (|A6p implies 

M 

A{u})fs{u}) = Ppa{uj)fa{uj). (A8) 

a=l 


Einally, /s(w) can be written as 

M , . 

Ol — 1 ^ ' 

M . . 

0 </.(-)<E^ (AlO) 

Oi — 1 ^ ^ 

0</s(w)<l, (All) 

where we used Eq. (IA6I1 and the fact that the Eermi-Dirac 
distributions satisfy 0 < fa(pj) < 1. The nonequilibrium 
distribution function is thus a linear combination of 
the Fermi-Dirac distributions of the reservoirs, and Eq. 
(lAllI) leads to an unambiguous definition of the local 
entropy for a nonequilibrium system given by Eq. m- 


Appendix B: Model of probe-sample coupling 


The scanning thermoelectric probe is modeled as an 
atomically sharp Au tip operating in the tunneling 
regime. The probe tunneling-width matrices may be 
described in general as = 2TrVnV^Pp{oj), where 

Ppipj) is the local density of states of the apex atom in 
the probe electrode and Vm,Vn are the tunneling matrix 
elements between the quasi-atomic wavefunctions of the 
apex atom in the electrode and the 7r-orbitals 

in the molecule. We consider the Au tip to be dom¬ 
inated by the s-orbital character and neglect all other 
contributions. The probe-system coupling is also treated 
within the broad-band approximation. The tunneling- 
width matrix describing the probe-system coupling 
is in general non-diagonal, and is calculated using the 
methods highlighted in Ref. 


Appendix C: Exact solution 


The exact solution to the Eq. (1311) with Ti —> 0 is 


Tp = T2 


2Ai 


where Ai is defined in Eq. (l40l) and simplifies to 


^ 7n^Tp^P{kBT2Y 

^ 20 Tpi 


a factor that appears in Eq. (155)) . 
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